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Abstract. The nearby accreting T Tauri star, TW Hya is viewed almost pole 
on, so its accretion and wind characteristics can be examined. A continuous 
set of echelle spectra of TW Hya, taken with the MIKE spectrograph on the 
Magellan2/Clay telescope at Las Campanas Observatory in April 2006 reveals 
systematic variations in the flux, velocity, and profile of the H-alpha emission 
line. These variations appear to be consistent with the photometric period of 2.8 
days. Absorption features recur at high outflow velocities (100 and 200 km s~^) 
in the wind. Additional spectra from 2004 show a similar repetitive pattern. 
This behavior suggests that: (a) accretion is not uniformly distributed over the 
stellar hemisphere in view; (b) stable structures are present in the chromosphere, 
most likely due to the stellar magnetic field configuration. Semi-empirical models 
of the atmosphere have been constructed to reproduce line profiles of H-a and 
He I, A10830 using the PANDORA code and to define the wind structure. These 
preliminary calculations suggest the mass loss rate is variable and comparable to 
H-a mass accretion rates in the literature, requiring a very efficient mechanism 
if the wind is powered only by accretion. 



1. Introduction 

The nearby accreting T Tauri star, TW Hya is oriented with its rotation axis 
almost along our line of sight, and the surrounding accretion disk approximately 
in the plane of the sky (Krist et al. 2000) providing a good opportunity to 
study its accretion and wind characteristics. The discovery of a hot, fast wind 
from this star (Dupree et al. 2005) calls for a determination of the mass loss 
rate, wind speed, and temperature in the outer atmosphere, and its relation to 
the mass accretion rate. A sufficiently robust stellar wind may lead to optical 
jets, could remove angular momentum from the star (Matt & Pudritz 2005), 
contribute to the opacity needed for X-ray absorption (Flaccomio et al. 2003), 
and influence the diminution of dust in accretion disks (Alexander et al. 2005). 
The chromospheric lines of H-a and He I A10830 are reported here for TW Hya 
to characterize the flux variations and line profiles. Detailed models of wind- 
sensitive line profiles of H-a and He I (A10830) are presented for expanding semi- 
empirical atmospheric models using the Avrett/Loeser PANDORA code for full 
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Figure 1 . Sample H-alpha profiles of TW Hya obtained at Magellan in April 
2006. The peak (Peak 1) used to measure the velocity shown in Fig. 2 {right 
panel) is marked. 
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Figure 2. Left: Flux in H-a as a function of the assumed photometric phase, 
Period=2.80 d (Lawson & Crause 2005). Zero phase is arbitrarily set at the 
beginning of the observations. The values are plotted twice. A sine function 
has been placed on the data to guide the eye. Right: Velocity of 'peak 1' in 
the H-a profile relative to the photosphere as a function of photometric phase. 
A sine function is overplotted. 



non-LTE effects. These models constrain both the atmospheric structure and 
the mass loss rate implied by the observed wind-scattered line profiles. 



2. Observations 

TW Hya was observed continuously for 3 successive nights in April 2006 at 
the Magellan/Clay telescope using the MIKE echelle spectrograph. Additional 
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Figure 3. Left: Gray scale representation of the Ha profiles of TW Hya for 
3 nights of observation in April 2006 where each spectrum has been divided 
by the average flux at each wavelength. The maximum Ha flux that occurred 
during night 2 is evident. Right: Gray scale representation of the Ha profiles of 
TW Hya where each spectrum has been divided by the maximum flux at each 
wavelength to display the absorption. Absorption at —100 km is replaced 
by a discontinuous jump to a new absorption feature at —200 km during 
Night 2 (~29 hrs). These absorptions are stable in velocity for hours, and not 
similar to the moving Discrete Absorption Components (DACs) observed in 
the winds of hot stars (cf. Howarth et al. 1995). 



optical spectra were taken previously in April 2004. The slit width of 0.75 
arcsec yielded a resolution of 36,000 near H-a. Infrared spectra of He I A10830 
were taken at KECK II using the NIRSPEC infrared spectrograph in May 2002 
and July 2005, and reported elsewhere (Dupree et al. 2005). Profiles of H-a 
showed substantial variation in flux and radial velocity over the course of all 
observations, (see Fig. 1, 2, and 3). The P Cygni profiles of He A10830 showed 
variations in the wind absorption among the 3 observations suggesting changes 
in the wind structure. The extent of absorption reaches ~ —300 km s~^ in both 
Ha and He I. Since H-a and A10830 are chromospheric lines, these velocities 
are supersonic and may be an indication of shocks and transient events. The 
photospheric escape velocity is ~500 km s~^, but at a distance of IR^ above 
the surface, the escape speed approaches 300 km s~^ so a small extension of the 
atmosphere could easily lead to mass loss. 
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Figure 4. Gaussian fits to the observed H-a profiles on the positive velocity 
side (left panel) reveal wind absorption for comparison to a model wind profile 
shown by the colored curve {right panel). Time dependent changes in the 
absorption indicate that multiple models are required with varying opacity 
and velocity structures. 
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Figure 5. Gaussian fits to the He, A10830 profile on the positive velocity 
side of a typical profile from 20 May 2002 {left panel) reveal wind absorption 
for comparison to the model wind profiles shown by the colored curve {right 
panel). As in the H-a line, the amount of the absorption varies, especially at 
low velocities, probably because of accretion. The stellar models for H-a and 
He I must be different which is not unexpected since the observations are not 
simultaneous. 
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Figure 6. Parameters for the Hydrogen Model (MODEL H-2004) and the 
Hehum Model (MODEL HE-2002). The temperature variation with radial 
height and the chromospheric acceleration are the principal differences be- 
tween the two models. The mass loss rate remains similar for both models. 
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3. Models 

The PANDORA code is a general purpose non-LTE atmospheric modehng and 
spectrum synthesis code (Avrett &: Locscr 2003). We invoke a model atmosphere 
which is spherical and expanding. PANDORA takes into account the time- 
independent optically-thick non-LTE transfer of line and continuum radiation 
for multilevel atoms and multiple stages of ionization including partial frequency 
redistribution. In calculating the line source functions we include the effects of 
atmospheric outflows. Multi-level atoms (3 levels for hydrogen, and 5 levels for 
helium, plus the continua) were used. 

The starting scmi-cmpirical models for the observations consist of temper- 
ature, density, and velocity and their radial dependence (see Fig. 6). In this 
calculation, the run of temperature and density are comparable to solar values; 
the density at 30,000K is consistent with values from the diagnostic C II lines 
(2326A) in TW Hya that indicate Ng ~ 4x10^ cm~^. The maximum coronal 
temperature is taken as the value (~ 10^-^ K) suggested by the emission measure 
distribution derived from CHANDRA X-ray spectra (Kastner et al. 2002). The 
velocity profile is constructed to match the velocities observed in the lines and is 
close to mass-conserving. The velocities and the radial height scales are varied 
in order to fit the observed absorption line profiles. 

The strengths of emission lines observed in T Tauri stars are substantially 
influenced by the hot plasma produced on the stellar surface as a result of 
the accretion flow. Thus we have fo cussed, not on reproducing the emission 
flux profile which is difficult to predict, but rather on the effects of the wind 
scattering on the line profile. 



4. Conclusions 

A continuous set of echelle spectra of TW Hya, taken with MIKE on Magellan2 
in April 2006 reveals systematic variations in the flux, velocity, and the profile 
of the H-a emission line. Absorption features recur at high velocity in the 
wind. Additional spectra from 2004 show a similar repetitive pattern. The 
flux variation indicates that accretion is not uniform and is consistent with the 
photometric variation found by Lawson and Crause (2005) which they attribute 
to the presence of accretion hot spots on the stellar surface. Absorption 'notches' 
appear on the short wavelength side of the H-a line at velocities near —100 and 
—200 km s~^. Similar notches appear in the He I A10830 profiles too. The 
fact that similar H-a profiles (with notches) occur in our data from 2004 and 
2006 indicates stable atmospheric structures. An obvious source creating such 
structures is a magnetic field that configures the chromosphere. Semi-empirical 
models of the atmosphere have been constructed to compute the scattering line 
profiles using the PANDORA code. The wind absorption appears to vary in both 
H-a and He A10830, and will require multiple models. The first approximations 
shown here for H-a and He suggest that the (spherical) mass outflow rate is on 
the order of 5 x 10""-'^'^ Mq yr"-*^. This value is comparable to the accretion rate 
indicated by H-a (Muzerolle et al. 2000) of 4x10"^'' Mq yr~^ and would appear 
to require a very efficient mechanism if the wind is powered only by accretion. 
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